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Abstract: We construct superfluid black hole solutions with two chemical potentials. 
By analogy with QCD, the two chemical potentials correspond to the baryon and isospin 
symmetries, respectively. We consider two systems: the back-reacted U{2) Einstein- Yang- 
Mills theory in 4-|-l dimensions and the 9-|-l-dimensional D3/D7 brane setup with two 
coincident D7- brane probes. In the D7-brane model, the identification of baryon and isospin 
chemical potential is explicit since the dual field theory is explicitly known. Studying the 
phase diagram, we find in both systems a quantum phase transition at a critical ratio of 
the two chemical potentials. However the quantum phase transition is different in the two 
systems: In the D3/D7 brane setup we always find a second order phase transition, while in 
the Einstein- Yang-Mills theory, depending on the strength of the back-reaction, we obtain 
a continuous or first order transition. We expect the continuous quantum phase transition 
to be BKT-like. We comment on the origin of this differing behavior in these apparently 
very similar models and compare to phenomenological systems. 
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1. Introduction 

In condensed matter physics, phenomenologically interesting systems are often strongly 
coupled. Famous examples are high Tc-superconductors and ultra-cold Fermi gases. Strong- 
ly correlated systems may undergo a phase transition at zero temperature which is driven by 
quantum fluctuations and thus named quantum phase transition [1] . Especially interesting 
are continuous quantum phase transitions which feature a quantum critical point. This 
quantum critical point influences the phase diagram also at non-zero temperature. In this 
influenced region, the quantum critical region, the system may be described by a critical 
theory even at finite temperature [2,3]. 

Gauge/gravity duality [4,5] provides a novel method for studying strongly correlated 
systems at finite temperature and densities. As such it should help to understand the 
systems described above [6-8] . Indeed remarkable progress was made in the application of 
gauge/gravity duality towards the description of superfiuids and superconductors following 
the results of [9], as well as of (non-) Fermi liquids [10]. 
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So far systems showing the transition to a holographic superfiuid have only been consid- 
ered with one control parameter, usually the ratio of chemical potential and temperature^. 
Naturally in such systems the phase transition is at a finite temperature and thus these 
systems have no quantum phase transition. In this paper we study a holographic system 
which shows a transition to a superfiuid phase and features an additional chemical poten- 
tial as a second control parameter. Since our starting point is a conformal field theory, we 
will construct two dimensionless control parameters: the ratio of one chemical potential 
to temperature and the ratio of the two chemical potentials. In such systems the ratio of 
temperature to one chemical potential at which the phase transition occurs may be tuned 
to zero by the second control parameter. Thus these systems may contain quantum phase 
transitions. 

Systems in which two chemical potentials can be tuned are often called imbalanced 
mixtures since two kinds of particles are present in imbalanced numbers. Examples are 
imbalanced Fermi mixtures where fermions with spin up and spin down are imbalanced [12], 
and QCD at finite baryon and isospin chemical potential where for instance up and down 
quarks are imbalanced [13] (see also [14]). Interestingly the phase diagrams of both these 
systems are very similar (see figure 1). In both systems there is a superfiuid state at 
low temperatures and at certain ratios of the two chemical potentials. In addition also 
the order of the phase transition agrees in both examples: At low temperatures (also at 
zero temperature) the transition is first order while at higher temperatures the transition 
becomes second order. Is it possible that there is an universal structure which relates these 
two different systems? 

In this paper we holographically study field theories which are expected to be similar to 
theories which describe the systems discussed above. Our field theories have a global U (2) 
symmetry which may be split into C/(l) x SU (2). This allows us to switch on two chemical 
potentials: one for the overall U{1) and one for a diagonal U{1) inside SU{2). In analogy 
to QCD, the chemical potential for the overall U{1) is the baryon chemical potential, while 
the one for the diagonal U{1) is identified with the isospin chemical potential. On the 
gravity side, we realize the U{2) gauge theory in two different ways: As a first model, we 
consider the U (2) Einstein- Yang-Mills (EYM) theory. In this model we allow the gauge 
fields to back-react on the geometry in order to get a coupling between the overall U{\) 
gauge fields and the SU{2) gauge fields. As a second model, we consider the D3/D7 brane 
setup with two coincident D7-brane probes which feature the U{2) gauge theory. In this 
model the interaction between the overall U{1) and the SU{2) gauge fields is obtained by 
the Dirac-Born-Infeld action. 

So far, holographic superfluidity in the Einstein- Yang-Mills theory in the presence of 
just an isospin chemical potential has been studied in e. g. [15-18]. In the probe approxima- 
tion [16], i.e. the gauge flelds do not influence the metric, a second order phase transition 
to a state is found which spontaneously breaks an Abelian symmetry. This spontaneous 
breaking creates a superfluid. In [18], the back-reaction of the gauge fleld on the metric 
has been added to this scenario. By increasing the back-reaction, the critical temperature 

^An exception is [11] where a second control parameter induced by multi-trace deformations is considered. 
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Figure 1: Phase diagrams of real world systems: (a) Imbalanced Fermi mixture in the canonical 
ensemble [12]: The spin polarization is the thermodynamic conjugated variable to the ratio of the 
chemical potentials favoring the different spins, (b) QCD at finite baryon and isospin chemical 
potential [13]. In both phase diagrams we observe a superfluid phase at small temperature and 
small ratio of two chemical potentials. In addition in both diagrams the phase transition is second 
order for large temperature and becomes first order at low temperatures. Both diagrams show a first 
order quantum phase transition. Both figures are reproduced by kind permission of the authors. 



decreases. Beyond a critical strength of the back-reaction, the phase transition is first 
order. There is a maximal value for the back-reaction beyond which the transition to the 
superfluid phase is not possible. 

The simple bulk action of the Einstein- Yang-Mills theory has the great virtue of being 
universal: The results may be true for many different dual field theories independently 
of their dynamics. Unfortunately this simple construction does not allow to identify the 
dual field theory explicitly. However it has been shown in [19-21] that the Einstein- Yang- 
Mills system can be embedded into string theory by considering the D3/D7 brane setup 
(see e.g. [22,23]). The dual field theory of the D3/D7 brane setup is known explicitly: 
It is = 4 Super- Yang-Mills theory coupled to M = 2 hypermultiplets. In this setup 
we work in the probe approximation, i. e. we consider N^'^ 1 D3-branes which generate 
the background metric AdS^ x and embed Nf = 2 D7-branes into the background 
space. The embedding of the D7-branes generates degrees of freedom which transform in 
the fundamental representation of the gauge group, the N = 2 hypermultiplets, which 
we denote as quarks in analogy to QCD. Here we have two quark flavors. Since the dual 
field theory is known explicitly, also our identification of the two chemical potentials as 
corresponding to the U{1) baryon and SU{2) isospin symmetries is explicitly realized. In 
this theory, mesonic bound states of the fundamental degrees of freedom are formed. The 
transition to the superfluid state is related to the condensation of vector mesons which 
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Figure 2: Sketch of the phase diagrams for the Einstein- Yang-Mills system for different strength 
of the back-reaction: In the white region the system is in the normal phase while in the blue region 
it is in the superfluid phase. The solid line marks a first oder phase transition and the dotted 
line a second order phase transition. In the normal phase at zero temperature the dual field theory 
contains an emergent one-dimensional CFT in the IR and the IR dimension of the operator depends 
on the ratio of the chemical potentials. For small back-reaction (a), the phase transition is second 
order for finite temperatures and we expect the quantum phase transition to be BKT-like. For 
intermediate back-reaction (b), there is a first order phase transition at large temperatures. At 
low temperatures the behavior is as for small back-reaction. For large back-reaction (c), the phase 
transition is always first order. Also the quantum phase transition is first order. 

spontaneously break an Abelian symmetry [19]. 

In both systems, i. e. in the EYM and in the D3/D7 model, the mechanism of breaking 
the Abelian symmetry is the same in the bulk: A non-zero vev of the time component 
of the gauge field At induces a chemical potential on the boundary theory. By fixing a 
gauge, we can choose the SU{2) gauge field in the direction of the third Pauli matrix to 
be non-zero, i. e. ^ 0. This breaks the SU{2) symmetry down to an Abelian symmetry 
which we call ?7(1)3. Beyond a critical value of the chemical potential, the systems become 
unstable against fluctuations of the gauge field pointing in some other direction inside the 
SU{2), for instance Al. (see e.g. [16] for the Einstein- Yang-Mills system and [24] for the 
D3/D7 brane setup). This instability is cured by the condensation of this gauge field A^. 
which then breaks the C^(l)3 symmetry. In the boundary theory the non-trivial profile of 
the gauge field A], induces a vev of the current (J^), but no source. Thus the breaking 
of the ^7(1)3 symmetry is spontaneous and the order parameter for the transition to the 
superfluid phase is given by (J^). For the D3/D7 brane setup we can explicitly write down 
the field content of the order parameter [19] , 

oc ipa^jxi' + (t>o-^dx<t) = ipulxipd + ipdlx-ipu + bosons , (1.1) 

where V' = ifpu,ipd) and cj) = i4'u,4'd) are the quarks and squarks duplet, respectively, a* 
denote the Pauli matrices and 7^ the Dirac matrices. 

In this paper, in addition to the time component of the gauge field inside the SU{2) 
A^, we switch on the time component of the U{1) gauge field At which induces the second 
control parameter, the baryon chemical potential. Tuning the two control parameters 
we can map out the phase diagram of both systems and find interesting similarities and 
differences (see figures 2 and 3 for a sketch of the phase diagrams) . In both cases the critical 
temperature where the phase transition occurs is finite at zero baryon chemical potential. 
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Figure 3: Sketch of the phase diagram for the D3/D7 brane setup: In the white region the system 
is in the normal phase while in the blue region it is in the superfluid phase. The critical temperature 
is first increasing and later monotonically decreasing as \ib I [i-i is increased. The phase transition 
is always second order. We also expect the quantum phase transition to be second order with 
mean-field exponents. There is also an emergent CFT in the IR at zero temperature in the normal 
phase. However, the IR dimension of the dual operator does not depend on ubI l^i and remains 
constant equal to 1. 

By increasing the baryon chemical potential, we can tune the critical temperature to zero 
and we obtain a quantum phase transition. However it is interesting that the details of 
the phase diagram are very different for the two systems, although they are expected to 
be dual to very similar field theories. For instance, the local as well as global symmetries 
match. The differences in the phase diagram are: In the Einstein- Yang-Mills theory (see 
figure 2) the critical temperature is monotonically decreasing as we increase the baryon 
chemical potential, while in the D3/D7 brane setup (see figure 3) the critical temperature 
first increases as the baryon chemical potential is increased. In addition in the Einstein- 
Yang-Mills setup, the system exhibits first and second order phase transitions depending 
on the strength of the back-reaction, while in the D3/D7 brane setup we obtain only second 
order phase transitions. Thus the question arises: What is the crucial difference between 
the systems which induces the different phase transitions? 

From the construction there is one obvious difference. In the Einstein- Yang-Mills 
system, the U{1) and SU{2) gauge fields only couple indirectly via the metric. In the field 
theory this means that the coupling of the currents which are dual to the gauge fields only 
occurs due to gluon loops. In the D3/D7 brane setup these loops are neglected due to 
the probe approximation. In this case the field theory currents directly interact with each 
other. These interactions are induced by the non-linear terms of the DBI action. Due 
to this difference it is understandable that the phase transitions may be different. The 
different couplings of the gauge fields to each other may lead to different RG flows and 
therefore to different IR physics which lead to differences in the phase diagram. 
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In addition we find an interesting difference in the origin of tlie quantum critical 
point in the systems. In the Einstein- Yang-Mills setup we can pinpoint the origin of the 
instability to the violation of the Breitenlohner-Freedman bound in an IR AdS2 region. 
This AdS2 region shows up as the near horizon region of the extremal Reissner-Nordstrom 
black hole. According to the AdS/CFT dictionary, the dual field theory thus contains a one- 
dimensional CFT in the IR (see figure 2). It is also important that the IR dimension of the 
dual operator depends on the ratio of the chemical potentials, such that the dimension can 
be tuned to an unstable value. In [25] it is argued that the violation of the Breitenlohner- 
Freedman bound will lead to a BKT-like phase transition. A common feature for this kind 
of transition seems to be the turning point in the phase diagram, such that the critical 
temperature slowly goes to zero as the ratio of baryon to isospin chemical potential is 
increased. In contrast to this behavior, the critical temperature in the D3/D7 brane setup 
goes to zero linearly. In this second model we do not obtain a violation of the Breitenlohner- 
Freedman bound in the IR, since the IR dimension of the dual operator does not depend on 
the ratio of the chemical potentials. Therefore we expect that the quantum phase transition 
is second order with mean field exponents. 

A similar difference occurs for the phase transition of chiral symmetry breaking via 
magnetic catalysis [26,27], if we compare this transition for the D3/D5 brane setup to the 
same transition for the D3/D7 brane setup [28,29]. In the D3/D5 model the quantum phase 
transition is BKT-like, while in the D3/D7 model it is second order. Here the difference 
between the two systems is more obvious: In the D3/D5 model both control parameters, 
the magnetic field and the baryon density, have mass dimension two. In the D3/D7 model 
the magnetic field still has mass dimension two but the baryon density has now mass 
dimension three such that the dimensions of the control parameters do not match. It 
is expected that the BKT-like transition only occurs if the two control parameters have 
the same dimension (see [30-33]). This is in contrast to field theories with two chemical 
potentials where the dimensions of the two control parameters match independently of the 
spacetime dimensions. This is a great advantage since these systems always satisfy the 
necessary condition for a BKT-like transition. However, this condition cannot be sufficient 
as our result in the D3/D7 brane setup shows. 

Comparing the phase diagrams obtained in our models (see figure 2 and 3) with the one 
obtained in imbalanced Fermi mixtures [12] and QCD at finite baryon and isospin chemical 
potential [13] (see figure 1), we see some similarities. In all cases the critical temperature is 
finite if the second control parameter, in our case the baryon chemical potential, is zero. By 
increasing the second control parameter we can tune the critical temperature to zero and 
we obtain a quantum phase transition. This seems to be a universal behavior for systems 
with two control parameters. However in imbalanced Fermi mixtures and QCD at finite 
baryon and isospin chemical potential shown in figure 1 the order of the phase transition is 
different from that in our models. In the models of figure 1, the phase transition is second 
order at large temperatures and becomes first order at low temperatures. On the other 
hand, in the holographic models this is different: For large back-reaction the behavior 
in the Einstein- Yang-Mills system is completely opposite. The phase transition is first 
order at large temperatures and becomes continuous at small temperatures. Also for small 
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back-reaction, we find a continuous quantum phase transition instead of a discontinuous 
one. 

Tlie difference in tlie order of the quantum phase transition may be related to the 
different behavior of the normal phase at zero temperature. For instance the BKT-like 
transition in the Einstein- Yang-Mills setup is possible since the theory is conformal in 
the IR and the IR dimension of the dual operator depends on ^s/fJ-i- In [25] a BKT- 
like transition has been discussed in conformal field theories. The transition may occur if 
two fixed points of the /3-function annihilate. The Einstein- Yang-Mills setup is the only 
one considered here which is conformal in the IR with tunable IR dimension of the dual 
operator, unlike both the models of figure 1 and the D3/D7 setup, such that a different 
order of the phase transition is plausible. In the D3/D7 probe brane setup we do not observe 
any change in the order of the phase transition which is always second order. Therefore, 
by comparing the different models, we conclude that the order of the phase transition is 
not universal and depends on the precise form of the interaction. 

The paper is arranged in the following way: In section 2 we study the back-reacted 
Einstein- Yang-Mills theory. In section 2.1 we present its action and equations of motion. 
In the normal phase the solution is given by the Reissner-Nordstrom black hole while in the 
superfluid phase the equations of motion are solved numerically by the shooting method. In 
section 2.2 we determine the thermodynamic quantities and construct the phase diagram 
in 2.3. In section 2.4 we discuss the phase diagram at zero temperature in detail and 
find an analytic expression for the quantum critical point. In section 2.5 we neglect the 
back-reaction of the Yang-Mills fields and analytically construct solutions in the superfluid 
phase for small baryon chemical potentials. This analytic solution determines the phase 
structure. 

In section 3 we investigate the D3/D7 brane setup. The embedding of the D7-branes 
is discussed in section 3.1. Its action and equations of motion are given in section 3.2. In 
section 3.3 we obtain the thermodynamic quantities and construct the phase diagram in 
section 3.3.1. Zero temperature solutions and the origin of the quantum phase transition 
are discussed in section 3.4. 

We conclude in section 4. 



2. Einstein- Yang-Mills Theory 



2.1 Action and equations of motion 

In this section we consider the U{2) Einstein- Yang-Mills theory in (4 -|- l)-dimensional 
asymptotically AdS space. The action is 



S 



d^xy 



-9 



1 

2k\ 



1 



1 



(2.1) 



^5 ^9mw ^9ym 

where is the five-dimensional gravitational constant, A = —12/R^ is the cosmological 
constant, with R being the AdS radius, ^mw the Maxwell and ^ym the Yang-Mills coupling. 
The U{2) gauge field is split into an SU{2) part with field strength tensor 

(2.2) 



rpa 
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where e is the total antisymmetric tensor and e = +1, and into an U{1) part with 
field strength tensor 



(2.3) 



The Einstein and Yang-Mills equations derived from the above action are 



where the Yang-Mills energy-momentum tensor Tf^i, is 



(2.4) 
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(2.5) 



Following [18], to construct charged black hole solutions with a vector hair we choose a 
gauge field ansatz 



A =<p{r)T^dt + w{r)T^dx , 
A =ip{r)dt. 



(2.6) 



The motivation for this ansatz is as follows: In the field theory we introduce a baryon and 
isospin chemical potential by the the boundary values of the time components of the gauge 
fields, (j) and ip. This breaks the U{2) symmetry down to a diagonal U{1) which is generated 
by r^. We denote this U{1) as U{l)s. In order to study the transition to the superfluid 
state, we allow solutions with non-zero (J^) such that we include the dual gauge field 
A]. = w in the gauge field ansatz. Since we consider only isotropic and time-independent 
solutions in the field theory, the gauge fields exclusively depend on the radial coordinate r. 
With this ansatz the Yang-Mills energy-momentum tensor in (2.5) is diagonal. Solutions 
with (J^) ^ also break the spatial rotational symmetry 50(3) down to S0{2) ^ such 
that our metric ansatz will respect only 50(2). Since the Yang-Mills energy-momentum 
tensor is diagonal, a diagonal metric is consistent. Following [18,34], our metric ansatz is 

ds^ = -N(r)a(rfdt'^ + -^^dr^ + r^f(r)-^dx^ + r^f(rf (dy^ + dz^) , (2.7) 

N[r) ' 

with N(r) = -2m{r)lr^ +r'^/R^. 

Inserting our ansatz into the Einstein and Yang-Mills equations leads to six equations 
of motion for m{r),a{r), f[r),(j){r),w{r),ip{r) and one constraint equation from the rr 

^Note that the finite temperature and chemical potential already break the Lorentz group down to 
S'0(3). 
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component of the Einstein equations. The dynamical equations may be written as 

" QNa^ ^ 6(72 V /' 6 ^ 

C = „ h cr — h 



N^a^ \r f N a 



(2.8) 



w(t? , fl 4/' A^' 



4^" = -i/j' 



3 _ a; 
r a 



(2.9) 



The equations of motion are invariant under five scahng transformations (invariant quan- 
tities are not shown), 

(/) a — )■ Ad, (j) ^ \(p , ip ^ Xip , 

{II) f^Xf, X-^w , 

(HI) r — > Xr, m, — > A^'m, w — t- Xw, (f> — > X(p, i/j — t- X^ , 

{IV) r Xr, m^X^m, R ^ XR, ^->A"V) V'-^-^^V) 

"YM — ^ AoyM , "MW — ^ AomW , 

(F) ^ Xip, auw A"^aMW , 

where in each case A is some real positive number. As in [18] we use (I) and (II) to set 
the boundary values of both a and / to one, so that the metric will be asymptotically 
AdS. Also we can use (III) to set r/j to one, but we will keep it as a bookkeeping device. 
We use (IV) to set the AdS radius R to one. The relation (V) allows us to set ccmw = 1 
by rescaling the baryon chemical potential, i. e. we can relate states with different baryon 
chemical potentials in different theories characterized by ayiw to each other. 

A known solution of the equations of motion is the AdS Reissner-Nordstrom black 
hole, 

= I^B-^ with Qi = mrl , 
= f{r) = l, (2.10) 

, 2(aYMgf + "MWgl) -.x, ^ «YMgf + «MWg| 

+ ^4 witi^ "^o = y + • 

In order to obtain the solutions in the superfiuid phase, i. e. w{r) ^ 0, we have to resort 
to numerics. We will solve the equations of motion using a shooting method. Wc will vary 
the values of functions near the horizon until wc find solutions with suitable values near 
the AdS boundary. We thus need the asymptotic forms of the solutions near the horizon 
r = r/i and near the boundary r — )• oo. 



4>{r) 


= IJ'i 




w{r) 


= 0, 


a{r) 


N{r) 




2mo 


= r2 
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Near the horizon, we expand ah fields in powers of et = r/rh — 1^1 with some 
constant coefficients. Three of these coefficients can be fixed as follows: We determine 
by the condition N{rfi) = which gives m{rh) = ^h/"^- Additionally, the time components 
of the gauge fields must be zero to obtain well-defined one-forms (see for example [35]), i. e. 
(j){rh) = and ilj{rfi) = 0. The equations of motion then impose relations among the other 
coefficients. A straightforward exercise shows that only five coefficients are independent, 



{alf^,w^,^1,i^l] , (2.11) 



where the subscript denotes the order of e^- All other near- horizon coefficients are deter- 
mined in terms of these five independent coefficients. 

Near the boundary, we expand all fields in powers of ef, = (rh/r)'^ <^ 1 with some 
constant coefficients. Again the equations of motion impose relations on these coefficients. 
There are seven independent coefficients 



{ml4,cp\,^Pl4,wl,f|] , (2.12) 



where here the subscript denotes the power of Cb- All other near-boundary coefficients 
are determined in terms of these seven independent coefficients. We used the scaling 
symmetries (2.9) to set o"q = /q = 1. Our solutions will also have Wq = since we do not 
want to source the operator in the dual field theory, i.e. the C/(l)3 symmetry will be 
spontaneously broken. In our shooting method we choose a value of cpQ = fij, the isospin 
chemical potential, and of -^q = ^s, the baryon chemical potential, and then vary the five 
independent near-horizon coefficients until we find a solution which produces the desired 
values at the boundary. 

In the following it will be often convenient to work with dimensionless coefficients by 
scaling out factors of r^. We thus define the dimensionless functions m(r) = m{r)/r^, 
(j){r) = (j){r)/rh, ■ilj{r) = ^{r)/rh and w[r) = w{r)/rh, while /(r) and a{r) are already 
dimensionless. 

2.2 Thermodynamics 

In this section we extract thermodynamic information from our solutions. The gravity 
solutions describe thermal equilibrium in the boundary field theory. In order to extract 
thermodynamic quantities from the gravity solutions we can use well-known methods of 
black hole thermodynamics. 

The temperature T in the boundary field theory is identified with the Hawking tem- 
perature of the black hole. The Hawking temperature for our black hole solutions is given 

by 

rjA A _ «YM ('^l)^ + «MW 

where k = sj d^j^^d^^ is the surface gravity of the black hole, with ^ being the norm of the 
timelike Killing vector, and in the second equality we write T in terms of the near-horizon 
coefficients. In the following we will often convert from the black hole radius to the 
temperature T by inverting the above equation. 



-] , (2.13) 
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The entropy S of the boundary field theory is identified with the Bekenstein-Hawking 
entropy of the black hole. For our ansatz we obtain 



S 



2^ . 



2Vnrl 



a 



YM 



12 K)^ 



(2.14) 



where Ah is the area of the horizon and V the spatial volume of the Minkowski space. 

The general statement of gauge/string duality which relates the field theory partition 
function to the string theory partition function may be used to calculate the thermody- 
namic potential of the boundary field theory, i. e. in our case the grand potential. In the 
gravity approximation, which we use in this paper, the grand potential Q is given as the 
temperature T times the on-shell bulk action in Euclidean signature. We thus analytically 
continue to Euclidean signature and compactify the time direction with period 1/T. We 
denote the Euclidean action as / and its on-shell value as /on-sheii- Since our solutions are 
always static, we can integrate out the time direction which produces an overall factor 
of 1/T. In order to simplify the expressions, we define / = I/T. From now on we refer 
to / as the action. I splits into three parts, a bulk term, a Gibbons-Hawking term and 
counterterms, 

I = Ihnlk + lGH + Ict- (2.15) 

The counterterms are needed to cancel the divergences of the bulk action and Gibbons- 
Hawking term which appear on-shell. To regulate these divergencies we introduce a hy- 
persurface at r = rbdy with some large but finite rbdy On the field theory side rbdy 
corresponds to an UV cutoff. Ultimately we will remove the cutoff by taking rbdy oo. 
Using the equations of motion, we obtain Ibuik'^^'^ ^"^^ ansatz 



fon-shell 
-'bulk 



V 1 



2 f2\/ 



rNa{r^f 



(2.16) 



=''bdy 



For our ansatz, the Euclidean Gibbons-Hawking term is 



fon-shell 
'GH 



' 3 



1 /• , V n f N' a- 

/ d^xVr V^n'^ = + - + 



(2.17) 



=rbdy 



where 7 is the induced metric on the r = rbdy hypersurface and n^dx^ = 1/ yj N{r) dr is 
the outward-pointing normal vector. The only divergence in the bulk action and Gibbons- 
Hawking term comes from the infinite volume of the asymptotically AdS space, hence, for 
our ansatz, the only nontrivial counterterm is 



Jon-shell = ^ / d^x^T = ^r^ViVa 



=rbdy 



Finally the grand potential is given by 



Vl = lim /on-shell • 

rbdy^OO 



(2.18) 



(2.19) 
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The baryon chemical potential is simply the boundary value of At{r) = ip{r) while the 
isospin chemical potential fij is the boundary value of ^{r) = (pir). The baryon charge 
density {Jt) and isospin charge density (Jf) of the dual field theory may be extracted from 
the on-shell action ion-sheii by 

(J) = - hm ^-^""-sheii ^ 27r^al,wT' <4m ('/'i ) ' + «MW (V'l ) \ \b 

(^l) ^ lb 

(2.20) 

Similarly, the current density (J^) is 

^ yrbd™ooMi(rbdy) Kl{a^f \ 12 {c^f ) 

(2.21) 

The expectation value of the energy-momentum-tensor of the dual field theory is [36, 37] 



^2 



/rr \ T ^ ^-^on-shell 

{Tij) = hm — .. = hm 



(2.22) 

-I ^•=ybdy 



where Z G {t,x,y,z} and = 1/2^/ N{r)dr^ij is the extrinsic curvature. We find 

For ^ = we recover the results obtained in [18]. Notice that the energy-momentum tensor 
is still diagonal such that the momentum is zero even in the superfluid phase where the 
current ( J^) is non-zero. This result is guaranteed by our ansatz for the gauge fields which 
implies a diagonal Yang-Mills energy-momentum tensor and a diagonal metric. 

For = 1/2 + {a\^,Ji] + g^,w/^l)/3, 4 = 4>i = 2/i/, = 2fiB, /| = 0, 
(j)^ = jli, and tp^ = jlB "we recover the correct thermodynamics for the Reissner-Nordstrom 
black hole, which preserves the 50(3) rotational symmetry. For instance, we find that 
{T^x) = (Tyy) = (T^^) and Q = -{Tyy). For solutions with non-zero (J^) the 50(3) 
symmetry is spontaneously broken to 50(2) and we find (T^x) 7^ (Tyy) = (Tzz)- However 
we also find 0, = — {Tyy) by just using the equations of motion as in [18]. 

Since the energy-momentum tensor is traccless (in Lorentzian signature), the dual field 
theory is scale invariant and describes a conformal fluid. The only physical parameters in 
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the dual field theory are thus the ratios hb/T, fJ-i/T and hb/ij-i- Since only two of them 
are independent from each other, we choose /U//T and fis/fJ-i to determine the physical 
state of the boundary field theory in what follows. 

2.3 Phase transition and phase diagram 

We expect a phase transition from the normal phase to a superfluid phase with a non-zero 
condensate (J^) as the baryon and isospin chemical potential are varied. From [18] we 
know that this phase transition occurs at zero baryon chemical potential. In the following 
we study the phase transition also at non-zero baryon chemical potential. 

Let us first map out the phase diagram of the U{2) EYM theory at finite temperature, 
baryon and isospin chemical potential for different values of the coupling constant oym- 
We start our discussion for small oym- Our numerical results are shown in fig. 4 and are 
confirmed by an analytic calculation at oym = presented in section 2.5. In the blue 
region the order parameter (J^) is non-zero and the system is in the superfluid phase while 
in the white region the order parameter (J^) is zero and the system is in the normal phase. 
We observe that the phase boundary moves monotonically to lower temperatures compared 
to the isospin chemical potential T/m as we increase the baryon chemical compared to the 
isospin chemical potential /^b////. The order of the phase transition does not depend on 
the baryon chemical potential and stays second order. At a critical value for the ratio of 
baryon to isospin chemical potential (/xb/^/)c we obtain a quantum critical point at zero 
temperature. In section 2.4 we determine this critical ratio analytically. Its value can be 
found in (2.29). 

By increasing ayu the area of the superfluid phase in the phase diagram decreases 
but the shape of the phase diagram stays the same until we reach a critical value for 
oym- Beyond the critical value (aYM)c,i = 0.365 ±0.001 we know from [18] that the phase 
transition to the superfluid phase becomes flrst order at zero baryon chemical potential. 
If we now increase the baryon chemical potential, we find a critical point where the phase 
transition becomes second order again (for a sketch see figure 2 (b)). The phase transition 
at zero temperature is still continuous and therefore a quantum critical point. If we increase 
ctYM) the critical point describing the change of the phase transition from first to second 
order moves to larger values of the ratio of baryon to isospin chemical potential. We find a 
critical value of ayyi where the zero temperature phase transition becomes first order and 
the quantum critical point disappears. Its value is given by (q!ym)c,2 = 0.492 it 0.008. For 
ttYM above this value the phase transition is always first order (for a sketch see figure 2 
(c)). 

2.4 Solutions at zero temperature 

In this section we consider the system exclusively at zero temperature. From the phase dia- 
grams presented above we see that for large baryon compared to isospin chemical potential 
the system is in the normal state. Since the normal state is described by a Reissner- 
Nordstrom black hole, the zero temperature limit is an extremal Reissner-Nordstrom black 
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0.05 




(a) (b) 

Figure 4: The phase diagram of the C/(2) Einstein Yang-Mills theory at finite temperature T, 
baryon /is and isospin chemical potential /xj for ayivi = 0.001 (a) and ayu — 0.1 (b): In the blue 
region the order parameter ( J^) is non-zero and the system is in the superfluid phase while in the 
white region the order parameter (J^) is zero and the system is in the normal phase. 



hole. Zero temperature is given by fixing the isospin chemical potential, 



^^I = . (2.24) 

As usual this extremal black hole features an AdS2 geometry in its near horizon region, 
i. e. in the IR. The full solution in the near-horizon region is given by 



ds2 = -12eW + ^+r2dx^ 



ayurh rh 



(2.25) 



where ^ = t — r^- According to the AdS/CFT dictionary, the dual field theory contains a 
one-dimensional CFT in the IR. Let us now consider this theory as we decrease the baryon 
chemical potential. From our numerical solutions we expect a phase transition towards a 
superfluid phase with non-zero vev (J^). This phase transition should be triggered by an 
instability in the normal state. In order to obtain this instability we consider fluctuations of 
the gauge field w{r) which is dual to the current about the extremal Reissner-Nordstrom 
background [38]. The equation of motion for this fluctuation is given by 



w" + -w'+ J^" Ji'^ w = ^, (2.26) 



where the prime denotes a derivative with respect to ^. This equation is the equation of 

36^2 



^ 3^2 ^2 ^2 

motion for a scalar field in AdS2 with effective negative mass squared mf^rr = L 2^ '2^ • 



Thus according to the AdS / CFT dictionary, the IR dimension of the dual operator can be 
tuned by changing the baryon chemical potential. Hence, the fluctuation is stable until the 
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mass is below the Breitenlohner-Freedman bound m^g = — 1/4.^ In our case, the bound is 
given by 

q 2 2 2 

< o • (2-27) 



^ayurh 2 

From this equation we may determine the baryon chemical potential at which the bound 
is saturated, 



t^B = . (2.28) 

With equation (2.24), we may determine the ratio between the baryon and the isospin 
chemical potential at which the Breitenlohner-Freedman bound is saturated. This ratio 
determines the point at which the system becomes unstable, 



^ * ^ . (2.29) 



Thus the Reissner-Nordstrom black hole may be unstable if oym < 1/ and a quantum 
critical point may exist if in addition omw is non-zero. This confirms our intuition ob- 
tained from our numerical results that at a given ratio of the baryon to isospin chemical 
potential a phase transition to a superfluid phase occurs. Unfortunately this calculation 
only determines the value for the ratio of baryon to isospin chemical potential where the 
system becomes unstable and not the phase boundary in general. For a continuous phase 
transition the two values coincide while for a first order phase transition, the transition 
always occurs before the instability is reached. Thus only for ayu < (oym)c,2 the phase 
boundary which is a quantum critical point and the critical value obtained here coincide. 
In [25] it is argued that the violation of the Breitenlohner-Freedman bound leads to a 
BKT-like transition. 

Naively we may assume that the superfluid phase is non-degenerate at zero temperature 
and the entropy is zero. In the gravity dual this is translated to a zero horizon radius of 
the black hole. The solution with zero horizon radius differs from the zero temperature 
solutions described by the extremal limit of the AdS Reissner-Nordstrom black hole with 
finite horizon size. Similarly to [17,39], we choose the following ansatz which is consistent 
with the numerical results near r = 0, namely 

(j) 4>i{r) , V'~V'i('")) u; ~ f^o + wi(r) , N r'^ + Ni{r) , 

m~mi(r), a ~ ctq + cJi(r) , /~/o + /i(r), 

such that all fields with index one go to zero at r = 0, e.g. /o + fi{r) — )■ /o as r — )• 0. 
Plugging the ansatz above in (2.8) and solving the equations of motion near the horizon 



^Note that at the boundary the geometry is ^4^5*5 where the Breitenlohner-Freedman bound is —4. 
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r = 0, we obtain the following solutions in the asymptotic forms 




2(3 

" r ^ N r i^!-^ — - — 5— , ~ Wo ' ^ 



^ = °' + ^j' ^~^°l^'-^2^^j ' 

(2.31) 

with 13 = JqWo. We can construct the full zero entropy solutions of the system by taking 
(2.31) as initial values near r = and integrate (2.8) numerically to the boundary using 
the shooting method. The result from that will describe the gravity dual of the superfluid 
ground state of the theory. 

It is important to note that a zero entropy solution is only consistent with ip being 
zero, i. e. no baryon chemical potential. Thus the domain walls we can construct from this 
asymptotics always have zero baryon chemical potential and coincide with the one found 
in [17]. At finite baryon chemical potential we expect a solution which interpolates between 
the domain wall solutions at zero baryon chemical potential and the extremal Reissner- 
Nordstrom solution in the normal phase. This solution should always contain a black hole 
with finite horizon radius and thus its entropy increases with the baryon chemical potential. 
By dimensional analysis we obtain S oc Vfj,^. 

2.5 The semi-probe limit 

In this section we study the EYM system taking just the back-reaction of the U{1) Maxwell 
field into account, i.e. oym = 0. We call this limit the semi-probe limit. From equation 
(2.29) we observe that there is the possibility of a quantum critical point at /is//-*/ = 
l/y/SaMW in this limit. In addition, the equations of motion (2.8) simplify significantly 
and we can obtain an analytical solution if we restrict ourselves to small baryon chemical 
potential fiB and small condensate (J^). The equations of motion in the semi-probe limit 
read ^ 



m 



(2.32) 



6 r 

since a = / = 1 if the back-reaction of the SU (2) Yang-Mills field is neglected. The 
equation for m and "ip can be integrated directly. 



m = "A + "MW/M f 1 _ 4 
2 3 V 

2 



/ ^2 



(2.33) 



Thus we are left with the two equations of motion for the SU{2) gauge fields in the given 
Reissner-Nordstrom background. By solving these equations numerically we can map out 
the phase diagram for ayivi = (see fig. 5). The phase diagram looks similar to the one 
where a small back-reaction of the SU{2) fields is included (see fig. 4). 
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Figure 5: Phase diagram in the semi-probe hmit: We compare the numerical data for the phase 
boundary (blue dots) with the analytic result (red line). We obtain a nice agreement for small 
baryon chemical potential where our approximation is valid. 

2.5.1 The expansion 

In the limit of small (ib and small ( J^), we can solve the equations of motion for and w 
analytically. For the case = 0, this has already been done in [40]. Similarly to [41], the 
solutions here are obtained as a double expansion in fiB and (J^) which are chosen to be 
proportional to the expansion parameters 6 and e, respectively. More precisely, we choose 
6 = fiB = ip{oo) and e = w\ (x (J^) from (2.21) where the tilde denotes dimensionless 
quantities which can be obtained by using (2.9) to set i? = r/j = 1. We make the following 
ansatz for (p and w 



Inserting the ansatz (2.34) into (2.32), we can construct a solution order by order in 6 and 
e. The possible solutions are restricted by certain boundary conditions. At the horizon 
r = 1, we demand that (p vanishes while w has to be regular. At the boundary, w is fixed 
to the expectation value (J^) tx e (2.21) while the isospin chemical potential jli associated 
to (j) receives finite corrections in 5 and e. 



w{r) 



Mr) + 5^ Mr) + 5^ M^) + ^ ('5') 
+ e' [(t>0,2{r) + 6^(l)2,2{r))+O[8^^) 
+ e^.^0,4(r) + O(<5V) 
+ O [e') , 

e (u;o,i(r) + 5^ W2,i{r) + 5^ Wi,i{r)) + O {d^e) 
+ e' {woM^) + 5^ W2,3{r)) + O {6^e^) 
+ eV5W + 0('5V) . 



(2.34) 
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The coefficient functions to lowest order, 

1 



„2 



(2.35) 



(l+r2)2 

are already known from [20] while the coefficients in the pure e expansion, that is with 
6 = 0, were first computed in [40]. To order e^, they read 

, , , / 1 \ 5 + 7r2 - 9r^ - 3r^ 71 

V r^J 96r^{l + r^) 6,720 

39 - 331r2 - 819r'^ - 369r6 + 156r2(l + r^f ln(l + 4^ ) 

= 20,160(1 +r2)5 ' ^'-''^ 

/ iVfl, / A 13 (-4,015,679 + 5,147,520 ln(2)) 
<^o,4(r) = /.,o,4 [^--2)+ ^o,4(r) , Wo,4 = 75,866,112,000 ' 

where $0,4 ('^) is a complicated function of r which we do not write down explicitly here. 
The IJ'im,n are determined by the regular boundary condition of w at the horizon r = 1 
and describe corrections to the critical isospin chemical potential /i/ = 4 at §"^e"' orders. 
The lowest order coefficient functions in the pure 5 expansion read 



(t>2,o{r) = fJ'12,0 (^i - , /i/2,0 = ^"MW (-17 + 24 ln(2)) 



Mr) = mfi (1 - ' (2.37) 
^^IA,o = ^ Omw [ - 5,495 + 864 vr^ ln(2) + 192 ln(2) (61 + 12 ln(2)2 - ln(8)) 
- 13,824Li3(l -i)- 13,824Li3(l + i) + 12C(3)] . 

For small baryon chemical potential jiB = 6, the critical isospin chemical potential for the 
phase transition will be corrected as 

= 4 + ^lnfl5^ + m,o'5' + 0(<J') • (2.38) 

This determines the phase boundary between the superfluid and the normal phase. We 
compare this analytic result with our numerical results in fig. 5. 
The lowest order coefficient functions in mixed orders read 

(l>2,2{r) = Ai/2,2 ( 1 - ^ ) + ^2,2{r) , 



/^/2,2 



(680,573 + 29,820 TT^ - 404,232 ln(2) - 1,406,160 ln(2)2) 



6,350,400 

2 / 13 + r2 (7 + 6r2 - 4 7r2(l + r^) + 24(1 + r^) \n{2f) 
W2,i{r) = "MW I 9(1^^2)3 (2.39) 

4 (3 - 20r2 + 2>r'^) ln(r) 2 (3 + 3r'^ + Ar"^ (-5 + ln(64))) ln(l + 



r2) 



+ 



9(l + r2)2 9(l + r2)2 

16r2 ln(r)2 8r2 (Li2(-r2) + Li2(l - r^) - 2Li2(l(l - r^)) ' 
3(l + r2)2 3(l + r2)2 



-18- 



where $2,2('") is a complicated function of r which we do not display explicitly here. 
2.5.2 The free energy 

Using the results from the last section, we compute the contribution to the free energy up 
to order (5"*e" for m + n < 4 from the gauge field term in the on-shell action 

a TpaAB 



A 



[ dr (r^{dr4>f - rN{r){drwf + -J—{4)wf] (2.40) 
Ji V N{r) J 

dr riV(r) {drw f j , 



25ym 

where V0I3 is the spatial volume of the field theory and /3 = l/T is the inverse temperature. 
In the q;ym — >■ limit, only ipir) contributes to the back-reaction which is described by 



YM -'1 

r^d) (d, 



N^r) = - 1 - il^-^a2,w<^' , (2-41) 



where the expansion parameter (5 = /xb is chosen to be small. 

For the background with vanishing condensate, i.e. a;(r) = 0, and 



^2 



- (4 + /X7o,2e^ + M/o,4e^ + + /X74,o5^ + M72,2(^^e^) , (2.42) 



the on-shell action is 

16 + 8/X7o,2e^ + (m7o,2 + 8/X7o,4) + 8fj,i2,oS'^ + {fJ'j2,o + 8a*/4,o) 
+ 2 im,2m,o + 4^/2,2) + O (<5f e^) 



_ pVohrl 

'-'vac — 2 

9ym 



(2.43) 



for p -|- g = 6. 

For the background where w ^0 has condensed, the on-shell action reads 

5sf = + 8M/o,2e' + (^/x2o_2 + 8^70,4 + ^]|^) + 8/^/2,0'^' 

+ (^12,0 + 8/X74,o) + 2 (/X70,2/X72,0 + 4^/2,2) Sh^ + O (,5^'e^) 
The difference in the values of the two on-shcll actions is 



(2.44) 



MP = S^-S.,= (-^e* + 0(^e.)) . (2.45) 

This result is known from [40] which remains robust in our back-reacted background. 

The free energy in the grand canonical ensemble is minus the value of the on-shell action 
times the temperature, hence the quantity AP determines the difference in the free energy 
between the normal and superfluid phase. In this case, the free energy of the superfluid 
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state is the smaller one because AP < 0, and this implies the stability of the superfluid 
phase. 

For small e and small 5, we have ~ (/2/ — jl'j{6))^. Using (III) in (2.9) to restore 
dimensions by taking \ = ^ T (2.13), the dimensionless /i/ will be replaced by ^ ~ ^ 
and thus we have ~ (Tc {5) — T). The cancelation of the term proportional to 5^e^ in the 
free energy difference suggests that the phase transition stays second order with mean field 
exponents as we increase the baryon chemical potential which coincides with our numerical 
result. 



3. D3/D7 Brane Setup 

3.1 Background and brane configuration 

In this section we investigate a string theory realization of the model studied above. We 
consider asymptotically AdS^ x spacetime which is the near-horizon geometry of a stack 
of D3-branes. The AdS^ x geometry is holographically dual to the = 4 Super Yang- 
Mills theory with gauge group SU{Nc). The dual description of a finite temperature field 
theory is an AdS black hole. We use the coordinates of [35] to write the AdS black hole 
background in Minkowski signature as 

d.^ = ^ (-Cdt' + /dx^) + (I)' {dg' + g'dnl) , (3.1) 

with dO^ the metric of the unit 5-sphere and 



o4 „4 



^4' ^y^' ^4 



f{g) = l-^, /(^) = 1 + ^, (3.2) 

where R is the AdS radius, with 

i?^ = AirgsNc a" = 2A a" . (3.3) 

The temperature of the black hole given by (3.1) may be determined by demanding regu- 
larity of the Euclidean section. It is given by 

T=^. (3.4) 

In the following we may use the dimensionless coordinate p = g/ gn^ which covers the range 
from the event horizon at /9 = 1 to the boundary of the AdS space at p — ?• oo. 

To include fundamental matter, we embed iVj coinciding D7-branes into the ten- 
dimensional spacetime. These D7-branes host flavor gauge fields A^ with gauge group 
U{Nf). This gauge field plays the same role as the gauge field in the Einstein- Yang-Mills 
systems. To write down the DBI action for the D7-branes, we introduce spherical coordi- 
nates {r, r^s} in the 4567-directions and polar coordinates {L,(j)} in the 89-directions [35]. 
The angle between these two spaces is denoted hy 9 {0 < 6 < vr/2). The six-dimensional 
space in the 456789-directions is given by 

dg'^ + g^dnl= dr^+r^dl^i+dL^ + L^dc/)^ = dg^ + g"^ {dO"^ + cos^ 0d(p'^ + sin^ 9dnl) , (3.5) 
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where r = gsinO, = r"^ + and L = gcosO. Due to the 5*0(4) rotational symmetry in 
the 4567 directions, the embedding of the D7-branes only depends on the radial coordinate 
p. Defining x = cos^, we parametrize the embedding by % = x{p) and choose = using 
the SO (2) symmetry in the 89-direction. The induced metric G on the D7-brane probes is 
then 



ds\G) 



dt^ + fdx^ + 



2R^\f ' J l-x^ 

The square root of the determinant of G is given by 



dg^ + R\l - X^)dnl (3.6) 



-G 



Q'ffil-x')\^l-x' + QHd,xy 



(3.7) 



where ^3 is the determinant of the 3-sphere metric. 

As in [24] we split the U{2) gauge symmetry on the D7-brane into U{1)b x SU{2)j 
where the U{1)b describes the baryon charges and SU{2)j isospin charges. As before we 
may introduce an isospin chemical potential m as well as a baryon chemical potential //^ 
by introducing non-vanishing time component of the non-Abelian background fields. Here 
we choose the generators of the SU{2)j gauge group to be the Pauli matrices cr* and the 
generator of the U{1)b gauge group to be = l2x2- This non-zero time-components of 
the gauge fields = At and break the U{2) gauge symmetry down to U (1)3 generated 
by the third Pauli matrix . In order to study the transition to the superfluid state we 
additionally allow the gauge field to be non-zero. To obtain an isotropic and time- 
independent configuration in the field theory, the gauge field A]j. only depends on p. This 
leads to a similar ansatz for the gauge field as in the Einstein- Yang-Mills theory, 



A = {At{g)(J^ + Al{g)(T^) dt + Al{g)a^dx . 



(3.8) 



With this ansatz, the field strength tensor on the branes has the following non-zero com- 
ponents, 

^ ox ^ XQ '^e X ' 



tx 



xt 



7 



(3.9) 



3.2 DBI action and equations of motion 

In this section we calculate the equations of motion which determine the profile of the 
D7-brane probes and of the gauge fields on these branes. The DBI action determines the 
shape of the brane embeddings, i. e. the scalar fields 0, as well as the configuration of the 
gauge fields A on these branes. We consider the case of Nf = 2 coincident D7-branes for 
which the non-Abelian DBI action reads [42] 



Sdbi = -Td7 Str J d^e VdetQ det (Pab [E^. + E^i{Q-^ - Sf^Ej,] + 2'Ka'Fab) 



(3.10) 
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with 



(3.11) 



and Pab the pullback to the Dp-brane, where for a Dp-brane in d dimensions we have 



,p, i, j = {p+l),...,{d- 1), E/,^ = 5^1, + Bf,t,. In our 



fi, = 0, . . . , {d — 1), a, b = 0, . 
case we set p = 7, d = 10, B = 0. As in [24] we can simphfy this action significantly by 
using the spatial and gauge symmetries present in our setup. The action becomes 



Sdbi = - Tm jd^^ Str^l det(G + 2to'F)| 



Td7 / d^e Str 



2 f^o\^ 



1 + G«G^^ (F3 )^ {a^ + 2F^^,J-,,aV0 + (J-,,)^ {a' 



2\2 / 2\2 



(3.12) 



where in the second line the determinant is calculated. Due to the symmetric trace, all 
commutators between the matrices a* vanish. It is known that the symmetrized trace 
prescription in the DBI action is only valid up to fourth order in a' [43,44]. However the 
corrections to the higher order terms are suppressed by NJ^ [45] (see also [46]). In [19,21] 
we used two different approaches to evaluate a non-Abelian DBI action similar to (3.12). 
First, we modified the symmetrized trace prescription by omitting the commutators of the 
generators cr* and then setting (cr*)^ = l2x2- This prescription makes the calculation of the 
full DBI action feasible. Second, we expanded the non-Abelian DBI action to fourth order in 
the field strength F. We obtained the same physical properties for the two approaches. We 
expect that the adapted symmetrized trace prescription also captures the relevant physics 
in this case such that we exclusively use the adapted symmetrized trace prescription in this 
paper. Using this prescription, the action becomes 



Td7 



d'? g'ffil - x') r^ip, X, A) + T2{p, X, A) 



(3.13) 



with 



'^i{p,X,A) 



i-x' + P\d,x? - 1(1 - x') {d,x,)' + l{i - x') (d,Al] ' 



P 



f 



7 



(l-x' + pHdpX?) (iX,-X2)Al 



(3.14) 



where the dimensionless quantities p = q/qh and A = {2T:a.')A/Qh are used. The fields 
Xi = At + and X2 = At — A^ are the gauge fields on the i-th brane. In [24] it is shown 
that the non-Abelian DBI action with A]. = decouples into two Abelian DBI actions 
in terms of these new gauge fields Xi. To obtain first order equations of motion for the 
gauge fields which are easier to solve numerically, we perform a Legendre transformation. 
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Similarly to [24, 35] we calculate the electric displacement pi and the magnetizing field 
which are given by the conjugate momenta of the gauge fields Xi and A^., 

In contrast to [24,35,47,48], the conjugate momenta are not constant any more but depend 
on the radial coordinate g due to the non-Abelian term A^A], in the DBI action as in [19,21] . 
For the dimensionless momenta pi and p]. defined as 



p ■■ 



P 



2Tra'TD7Qi 



(3.16) 



we get 



Pi ~ ^r^ / In 1 Px ~ 



2fTi{p,x,A) ' 2 
Finally, the Legendre-transformed action is given by 

Sdbi = 5'dbi — / d*.^ 



+ 



Ti{p,x,A) T2{p,x,A) 



5 



J d«e Q^fRl - X^)^^-X^ + p'{dpX? V{p, X, A,p) , 



(3.17) 



(3.18) 



with 



V{p, X, i, P) = ( 1 - ^^^^ {{Xi - X2)Al) 



'1 + 



8(pi)2 



+ Jl + 



8(P2)2 



8(pi)' 



pW(i-x2)3 

(3.19) 



This action agrees with the one for finite baryon and isospin chemical potential (see [24]) 
after — t- and with the one for the superconducting state at pure isospin chemical 
potential (see [19]) after pi — > —p2 and p^ ^fPx- The change in p}^ has to be done such 
that the definitions agree in both cases. 

Then the first order equations of motion for the gauge fields and their conjugate mo- 
menta are 



dpXi 



dpAl 



^-^^^^^^^^§^m{p, X, ApniP, X, A,p) , 

p-^f^i^-x^r 

2^i-x' + pHd,xY . 

-pJV{p,x,A,p) , 



p'fil - x2)2 

dpPi/2 = ± 

/(I - x')V^-x' + pHdpXyi' 



/(I - x')Vl -X' + P^[dpx)W / ,^ 
87r'^pfW{p,x,A,p) \ ^) ^ ' 



(3.20) 



dpPx = 



87r^pfW{p,x,A,p) 
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with 



Ui{p,x,Ap) 



w{p,x,Ap) 



1 + 



1 



p6/3(l_^2)3 



(3.21) 



8(pi)^ 

p6/3(l_^2)3 



1 + 



+ W1 + 



p6/3(l_^2)3 



For the embedding function % we get the second order equation of motion 



do 



p'ffil-X^)idpx)V 



p^ffx 



v/i-x' + p'(^px)U 

24{l-x'+pHd,xr) 



[3{l-x')+2p\d,xf]V 



p'fHi-x 



2^3 



W 



(3.22) 



We solve the equations of motion numerically by integrating them from the horizon at 
p = 1 to the boundary p = oo. The initial conditions may be determined by the asymptotic 
expansion of the gravity fields near the horizon 



X; 



Pl/2 = Ol/2 



pI 



X = Xo 



bi 



327r2 \ Bi B2 



{Bi + B2) {p - 1 
3X0 / 



If 


+0{{p- 






+0{{p- 


If), 


I? 


+0{{p- 


If), 




+0{{p- 


1)^), 


If 


+0{{p- 


1)^), 



(3.23) 



with 

The terms in the asymptotic expansions are arranged according to their order in (p — 1). 
There are four independent parameters {a,fci,&25Xo} which have to be determined. In 
order to obtain the field theory quantities we determine the asymptotic expansion of the 



(3.24) 
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gravity fields near the AdS boundary 



Pi 



X 



di 



= l^i 

di 
-dl 



+ 



m 
P 



(3.25) 



Note that the factor of two in is consistent with the earlier definitions in [21] since here 
we have a different definition of the conjugate momenta (factor Nf). In this asymptotic 
expansion we find seven independent parameters {/xj, di,d\., m, c}. Using the transformation 
of the gauge field from to At and A^, the independent parameters of the gauge fields 
Xi, {fii, di}, may be translated into parameters of the asymptotic expansion of At and , 



A/ -- 



A? 



^(Ai + A2) 

:(Ai - A2) 



dB = df = dl + d2 , 
di = df = dl — d2 ■ 



(3.26) 



These parameters may be translated into field theory quantities according to the AdS/CFT 
dictionary (for details see [35]): /xs is the baryon chemical potential, pi the isospin chemical 
potential, 



P-B 



2 PB 

A T ' 



PI 



2/X7 

A T ' 



(3.27) 



the parameters d are related to the vev of the flavor current J by 



dB = d\ 



di = dl 



dl 



and m and c to the bare quark mass Mq and the quark condensate (V'V')) 



(3.28) 



m = 



2Mq 



VXNfNcT^ ' 



(3.29) 



respectively. There are three independent physical parameters, e.g. m, ps and pi in 
the grand canonical ensemble. The asymptotic expansion close to the horizon has four 
independent solutions. These parameters may be fixed by choosing the three independent 

physical parameters, i. e. the state in the field theory and by the constraint that Al. goes to 
zero at the boundary, i. e. the ^7(1)3 symmetry is spontaneously broken. We use a standard 
shooting method to determine the parameters at the horizon. 
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3.3 Thermodynamics 

In this section we study the contribution of the D7-branes to the thermodynamics. Ac- 
cording to the AdS/CFT dictionary the partition function Z of the boundary field theory 
is given in terms of the Euchdean on-shell supergravity action /on-sheii> 

Z = e-^°"-hcii . (3.30) 

Thus the thermodynamic potential, i. e. the grand potential in the grand canonical ensem- 
ble, is proportional to the Euclidean on-shell action 

Q = -TlnZ = r/on-shell ■ (3.31) 

To calculate the contribution of the D7-branes to the grand potential, we have to determine 
the Euclidean version of the DBI-action (3.13) on-shell. For this purpose, we first perform a 
Wick rotation in the time direction. Next we renormalize the action by adding appropriate 
counterterms Id (see [49] for a review and [50] for probe D-branes). In our case the 
counterterms are the same as in [24,35,51], 



XN.NfVsT'^ 

-'rt — 

128 



Pmax-"^^)^ -4mc , (3.32) 



where pmax is the UV-cutoff and V3 the Minkowski space volume. Then the renormalized 
Euclidean on-shell action may simply be written as 



^R = - 



32 \ Nf 4 



(Pmin -in'^)'^ - Amc ^ , (3.33) 



where pmin determines the minimal value of the coordinate p on the D7-branes, i. e. pmm = 1 
for black hole embeddings which we consider exclusively in this paper and 

G{m,il)= r dp[p'ff{l- x") (jiip,x,A) + T2{p,x,A)) - Nf {p^ - pm)] . (3.34) 

Pmin 

In the following we consider the dimensionless grand potential W7 defined as 

= (3.35, 

By considering the variation of the grand potential with respect to the gravity fields, it can 
be shown (see [21, section 5.3] for the pure isospin case) that the above definition of the 
grand potential is consistent and that the order parameter d], is not a thermodynamical 
variable. 



3.3.1 Phase transition and phase diagram 

We expect that a phase transition occurs between a normal fiuid phase and a superfiuid 
phase. At zero baryon chemical potential we know from [19,21] that the phase transition 
from the normal phase to the superfiuid phase is second order with mean field exponents. 
In this section we consider the phase transition at non-zero baryon chemical potentials. 
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First we map out the phase diagram of the given theory with zero quark mass m = 
0. The phase diagram is shown in fig. 6. As we increase the baryon chemical potential 
the transition temperature to the superfluid phase first increases. For /x_b/a*/ ^ 0-4 the 
transition temperature monotonically decreases to zero as the baryon chemical potential 
grows. We can show numerically that the phase transition is always second order. In 
the next subsection we will show numerically that the transition temperature is zero at 
fJ-s/fJ-i ^ 1.23, and hence we obtain a quantum critical point. 




0.0 0.2 0.4 0.6 0.8 1.0 1.2 

fJ-B 



Figure 6: The phase diagram for fundamental matter in thermal strongly-coupled J\f = 2 SYM 
theory at zero quark mass with /x/ the isospin chemical potential, fiB the baryon chemical potential, 
T the temperature and A the 't Hooft coupling: In the blue region the order parameter ( J^) is non- 
zero and the system is in the superfluid phase while in the white region the order parameter (J^) 
is zero and the system is in the normal phase. 



3.4 Zero temperature solutions and quantum critical point 



In this section we study the D3/D7-brane setup at zero temperature. The zero temperature 
limit is given by p/j — )• 0, i. e. / = / = 1. The induced metric on the D7-branes may now 
be written in {L,r) coordinates (see equation (3.5)) 



+ ^ ^-dt^ + d.-)+^(l + (5.L)^)dr^ + ^...3 



The square root of the metric is now 



-G 



+ {drLY , 



(3.36) 



(3.37) 



and using the adapted symmetrized trace prescription, the DBI action becomes 

^3 



'S'dbi = —Td7 / d^ 



i{f,L,A) + E2if,L,A) 



(3.38) 
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with 



1\2 



1 + {dfLY - 2{dfXiY + 2{dfAi,) 



7 



27r2(f2+L2)2 

where the dimensionless quantities are now defined by 



1 + {dfif) f (Xi - X2)Al 



(3.39) 



r ~ L ~ 2-Ka . 

r = — , L = —, A = A. 

R R R 



(3.40) 



In the normal phase, i. e. A}^ = 0, the equations of motion for the gauge fields Xi for the 
massless embedding L = 0, 

dfXi = ' (3.41) 

'f6 + 8d2 

can be solved analytically [47]. The solution expressed in terms of incomplete Beta func- 
tions is given by 

(]^ I 11 



Prom the asymptotic form near the boundary, we can read off the chemical potential and 
the density (see equation (3.25)), 



'~ Til) 



(3.43) 



In the normal fluid phase we may now consider fluctuations Z± = A^. it 1^2 (see [24, 52] 
for more details) and look for an instability which may lead to a phase transition. For the 
flat embedding L = 0, the equation of motion for the fluctuation Z± at zero momentum is 
given by 



d'^Z+{r) F'{f)dZ+{r) 



df2 ~^ Fif") df 
d2Z_(f) F'(r)dZ^(r) 



7 



df2 



F{f) df 



A^/2-K 

7 

4^/27r 



^+(f) = 0, (3.44) 



X1-X2 Z^{f) = 0, (3.45) 



with Co = \l ^ a'uj coming from the ansatz A\.{r, t) = Al^{f)e The analytical expression 



for Xi{f) is given in (3.42) and 

/ 

F{r) = 

V 



1 - 2 (dfXi 



+ 



l-2(dfX2 



(3.46) 



The system is unstable if the imaginary part of the quasinormal frequency is positive. For 
a massless embedding L = 0, we find this instability at /is//// ~ 1.23. 
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3.4.1 What is the origin of the instability? 

In the back-reacted Einstein- Yang-Mills theory we see that the instability of the extremal 
Reissner-Nordstrom black hole which triggers the phase transition to the superfluid phase 
can be related to the violation of the Breitenlohner-Freedman bound of the field A], in 
the near horizon geometry of the black hole. This near horizon geometry features an 
AdS2 factor (see section 2.4). In this section we investigate the instability in the D3/D7 
brane setup at zero temperature. In the following we show that in contrast to the back- 
reacted Einstein- Yang-Mills system, the instability in the brane setup cannot be related 
to the violation of the Breitenlohner-Freedman bound in the far IR, i. e. the fiuctuation 
of the field A^. does not become unstable in the far IR. In this system we expect that the 
instability occurs in the bulk. 

Let us now consider a fluctuation Al. about the solution in the normal phase given by 
equation (3.42) for zero quark mass, i.e. L = 0. The equation of motion for this fluctuation 
is given by 

(Al)"^^(AlY + f%^Ai = 0. (3.47) 

As in section 2.4 we consider the equation of motion in the far IR. Since there is no horizon 
at zero temperature in the brane setup, the expansion is around f = 0. The equation of 
motion becomes trivial, d-Al. = 0, if we expand to flrst order only. Thus Al./r satisfies the 
equation of motion of a massless scalar in AdS2. In contrast to the EYM setup, the mass 
cannot be tuned by changing a UV quantity and the Breitenlohner-Freedman bound can 
never be violated^. Hence the fiuctuation is stable in the IR. We expect that the instability 
observed by the numerical study of the quasinormal modes above occurs along the flow. 
Thus the origin of the instability in this brane setup is different compared to the one in 
the back-reacted Einstein- Yang-Mills system studied in section 2.4. 

4. Conclusion 

We have considered holographic models of field theories with global U (2) symmetry. The 
C/(2) symmetry allows us to switch on two chemical potentials: a baryon and an isospin 
chemical potential. Holographically we realized the global U{2) symmetry in two different 
ways: first by an U{2) Einstein- Yang-Mills theory and second by the D3/D7 brane setup 
with two coincident D7-branes. We mapped out the phase diagrams for both systems and 
found interesting similarities and differences which we already discussed in the introduction. 

It would be interesting to study the origin of these differences in the phase diagrams in 
more detail. For example a detailed analysis of how the order of the phase transition can be 
changed by varying the form of the interaction would be attractive. In addition the study of 
back-reaction effects in the D3/D7 model may lead to new behavior in the phase diagrams. 
Furthermore it is important to establish a full understanding of the instability mechanism 
in the D3/D7 brane setup and its difference to the violation of the Breitenlohner-Freedman 
bound found in the Einstein- Yang-Mills theory. This may lead to a characterization of the 
universality classes of quantum phase transitions. 
*We thank K. Jensen for pointing this out. 
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